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We compute the absorption spectrum of strongly repulsive one-dimensional bosons in a disordered 
or quasi-periodic optical lattice. At commensurate filling, the particle-hole resonances of the Mott 
insulator are broadened as the disorder strength is increased. In the non-commensurate case, map- 
ping the problem to the Anderson model allows us to study the Bose-glass phase. Surprisingly we 
find that a perturbative treatment in both cases, weak and strong disorder, gives a good description 
at all frequencies. In particular we find that the infrared absorption rate in the thermodynamic limit 
is quadratic in frequency. This result is unexpected, since for other quantities like the conductivity 
in one dimensional systems, perturbation theory is only applicable at high frequencies. We discuss 
applications to recent experiments on optical lattice systems, and in particular the effect of the 
harmonic trap. 
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I. INTRODUCTION 

In recent years, developments in the field of ultracold 
atomic gases have considerably enlarged the possibili- 
ties for exploring the physics of strongly correlated sys- 
tems pp. For instance, the study of quantum phase tran- 
sitions, a subject of continuous theoretical interest, has 
been strongly stimulated by the experimental observa- 
tion of the superfluid to Mott insulator transition using 
optical lattices [2]. Indeed, ultracold atom systems offer 
us an unprecedented control over the system parameters, 
which can allow us to ultimately understand the physics 
of very hard problems such as the phase diagram of the 
Hubbard model in two-dimensions [3]. 

A particularly interesting and fertile arena is the study 
of disordered ultra-cold atoms. Random potentials can 
be introduced in a controlled way by laser beams gener- 
ating speckle patterns 01 El El [Zl IH1 E] , or by loading in 
an optical lattice a mixture of two kinds of atoms, one 
heavy and one light. When the heavy atoms become ran- 
domly localized in the lattice, they will act as impurities 
for the lighter atoms [TU]. Another available technique 
is to superimpose two optical lattices with incommen- 
surate periodicities, thus generating a quasi-periodic po- 
tential [TT]. The quasi-periodic lattices as well as the 
speckle patterns have been recently used in the experi- 
mental efforts to observe the effects of Anderson localiza- 
tion in dilute Bose gases expanding in highly elongated 
traps [T21 [15] , Since Anderson localization is a single- 
particle effect, the next logical step is to study the inter- 
play of disorder and interactions in strongly interacting 
ultra-cold atomic systems [131 UHl HS1 HZ]- The latter 



may be accessible by tuning intcrparticle interactions us- 
ing Feshbach resonances [TH] , or by loading the atoms in 
sufficiently deep optical lattices [5] and/or strongly con- 
fining them to low-dimensions in tight traps pp. 

In the context of the efforts described above, one of the 
experimental challenges is to observe in ultra-cold gases 
clear signatures of the theoretically-predicted transition 
from a superfluid to the Bose-glass phase [El [20j 121"] . 
A pioneering step in this direction was recently taken 
by the Florence group by using lattice modulation spec- 
troscopy This technique consists in heating an ultra- 
cold gas loaded in a optical lattice by periodically mod- 
ulating the depth of the lattice [22J H3 ij. When per- 
turbed in this way, the gas is driven out of equilibrium 
and absorbs energy. When the perturbation is switched 
off, and after re-thermalization, the broadening of the 
momentum distribution around zero momentum is taken 
as a measure of the energy absorbed by the system dur- 
ing the lattice modulation [22 . On the theory side, 
for non-disordered lattices, the calculation of the energy 
absorption rate due to the lattice modulation was first 
performed analytically within linear response theory by 
some of the present authors [23 . These results were con- 
firmed and extended beyond linear response using time- 
dependent density-matrix renormalization group meth- 
ods, both for the case of bosons [21] and fermions [25] . 
More recently, for disordered optical lattices, the energy 
absorption rate has been numerically calculated using full 
diagonalization in small systems [251 [27]. Other meth- 
ods that have also been discussed in the literature for 
detecting signatures of the effect of disorder (or quasi- 
periodicity) on interacting boson systems in one dimen- 
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sion focus on the momentum distribution [TSJ [55] in dis- 
ordered systems as well as quasi-periodic systems or the 
expansion dynamics in quasi-periodic systems [15] . 

Let us consider a one dimensional disordered Bose gas 
described by the following Hamiltonian: 



H = -J£(&j +1 fe, + h.c.) + n ^ !) 



(i) 



where bj denotes the boson annihilation operator at sites 
j t rij — b^bj being the local density. Here J and U are 
the usual parameters of the Bose-Hubbard model corre- 
sponding to the tunneling rate and the onsite repulsion 
(U > 0). The last term in the rhs of Eq. (HI accounts 
for the presence of both the harmonic trap, Vj°, and the 
disorder potential tj. 

The distribution of the on-site energies €j in Eq. (lj) 
depends on the specific choice of the random (or pseudo- 
random) potential. In this work we extensively compare 
two cases: 

a) ej uniformly distributed in [—A, A] (2) 

b) ej — A cos(27rjer), a irrational, (3) 

where A measures the strength of the disorder. From 
the experimental point of view, speckle patterns and lat- 
tice containing heavy atom impurities can be modeled by 
case a) whereas quasi-periodic potentials obtained by su- 
perimposing two optical potentials with incommensurate 
periodicities d\jdi = a are described by case b). 

In this work we assume that the hopping amplitude 
J in Eq. ([lj is modulated periodically in time according 
to J(t) = J + <5 J cos u>t, where u> is the modulation fre- 
quency. We calculate the energy absorption rate within 
linear response theory, which is valid for weak lattice 
modulations 6 J <C J. 

We shall restrict ourselves to the strongly repulsive 
regime of Eq. ([lj, where the on-site repulsion is large 
compared to both the hopping amplitude and the dis- 
order strength (U J, A). We first discuss systems in 
the thermodynamic limit by setting Vj 10 = 0. Effects of 
the harmonic trapping potential will be discussed later, 
in Sect|V[ We also set H = 1 to simplify the notation. 

In the thermodynamic limit both the number N of 
bosons and the length M of the chain diverge. The ratio 
v = N/M is instead finite and corresponds to the filling 
factor, i.e. the average number of bosons per lattice site. 
We must distinguish between two physically different sit- 
uations. For incommensurate fillings v < 1, the system 
has gapless excitations. Thus, a good approximation to 
the absorption rate at frequencies lo <C U can be obtained 
by formally taking U — > +oo and mapping the result- 
ing hard-core bosons to non-interacting fermions using a 
Jordan- Wigner transformation (see e.g. Ref. |29j. The 
resulting single particle problem can be easily solved for 
a given choice of ej and V ho and the energy absorption 
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FIG. 1: Sketch of the full absorption spectrum at incom- 
mensurate filling. Here W ~ max(4 J, 2A) is the effective 
bandwidth. The peak at to ~ U stems from particle-hole ex- 
citations, whereas the low frequency absorption appears as a 
consequence of the formation of a Bose-glass at incommensu- 
rate filling. 



obtained. On the other hand, for unit filling (y = 1) the 
ground state is a Mott insulator with a gap. In this case, 
the elementary excitations are particles and holes [23] 
and in a homogeneous system the absorption can only 
occur at frequencies w ~ U . In the absence of disor- 
der, the response of the system to the lattice modulation 
has been computed using degenerate perturbation theory 
within the subspace of particle and hole excitations [23 . 
The same methods can be generalized to deal with dis- 
ordered CctSC . MS shown below. 

The resulting general picture of the energy absorption 
is depicted in Fig. [T] The low frequency < u> < W, 
where W ~ max(4 J, 2A) is the effective bandwidth, is 
only present for incommensurate fillings, v < 1. A second 
distribution is centered at frequency u — U and comes 
from particle-hole excitations generating doubly occupied 
sites. The width of this absorption line is also given by 
W. 

The paper is organized as follows. In Section [n] we 
present the general formalism needed to calculate the 
absorption rate for incommensurate and unit filling. In 
Section [III] and Section [TV] we present our analytical and 
numerical results obtained for the random and the quasi- 
periodic potential, respectively. In Section [V] we discuss 
the effects of a trapping potential. Finally in Section |VI| 
we provide our conclusions. A derivation of the general 
formula [see Eq. d9|] for the absorption rate valid for weak 



disorder is given in the Appendix VII ) 



II. 



ENERGY ABSORPTION: LINEAR 
RESPONSE THEORY 



For weak perturbations, corresponding to 6 J <C J, the 
energy absorption rate can be calculated using linear 
response theory. The general formula has been first de- 
rived in Ref. [25] and in the presence of disorder it takes 
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the form 



E^ = -5J wlm[-xif(w)], 



(4) 



where is the Fourier transform of the retarded cor- 

relation function x*r(*) = - i0 (*)([-^(*)> -^(°)]) of the 
hopping operator K = +h- c -)' being 6(i) 

the step function. In Eq . (E| the bar means average over 
different disorder (Sect. IV I or quasi-periodic (Sect. Ill) 
realizations. 

In general, the calculation of the correlation function 
in Eq. Q is a complicated many-body problem. How- 
ever, in the limit of strong repulsion where U 3> J, A, 
calculations are considerably simplified by the fact that 
we can accurately restrict ourselves to work within a sub- 
space of the total Hilbcrt space, whose detailed structure 
depends on the filling. In the case of an incommensu- 
rate filling (i.e. v < 1) this subspace can be described 
in terms of non-interacting fermion states (that is, Slater 
determinants). For commensurate filling (i.e. unit fill- 
ing, v = 1) we can restrict ourselves to the subspace with 
one particle and one hole excitation, as described below. 



Incommensurate filling 



In the previous expressions e a and ip a (j) are the eigen- 
values and eigenfunctions of H', respectively. In Eq. ^ 
/(e) = (exp[(e — /i)/T] + l) _1 is the Fermi-Dirac distribu- 
tion function at a temperature T and chemical potential 
[i. The latter is fixed by the normalization condition 
v = f(e a ). At zero temperature the only relevant 
processes in Eq. ([7]) correspond to transitions from an 
occupied level (with energy e a < (x(T — 0)) to an un- 
occupied level (sp > /jl(T = 0)). In particular, for unit 
filling the absorption ^ vanishes, consistently with the 
fact that a Mott insulator can only absorb at much higher 
frequencies u> ~ U. 

In the absence of disorder [i.e. for e, = in Eq. dfjj], 
the hopping modulation commutes with the Hamiltonian 
H' and therefore the absorption rate vanishes to all or- 
ders, even beyond linear response. In the above expres- 
sion, this is reflected in that, for ej = 0, the eigenstates 
of H' become plane waves, ipk(j) °c e J , with energy dis- 
persion Ek — — 2 J cos k (k being the lattice momentum). 
Thus the matrix Q is diagonal, i.e. tCkk' — ^k,y c °s 2 
which, together with the factor /(£&) — f(sk') in Eq. ^ 
makes E u vanish. 

At weak disorder (i.e. J> A), the absorption rate ^ 
can be evaluated using perturbation theory (the details 
can be found in the Appendix), which yields: 



For filling v < 1 and large on-site repulsion J,A<£lU, 
we take the hard core limit U — > +oo. Bosons are then 
mapped onto non interacting spinless fermions via the 
Jordan- Wigner transformation: 



exp 



' 3-1 
. fe=l 



(5) 



where 



satisfy fcrmionic commutation relations 



{cpCj} — and {c_-,ct} = 1. Under the above transfor- 
mation, the Hamiltonian ([l]) is mapped onto the single 
particle Hamiltonian: 



^ = - J E(4 + i c ,+ h - c -) + E e 3 c ]< 



(6) 



where the on-site repulsion U has disappeared and the 
hopping operator K in Eq. Q is now given by K' — 
— J2j(Cj +1 Cj + h.c). Notice that the mapping of observ- 
ables that are non local in space is far less trivial, an ex- 
ample is the momentum distribution studied in Ref. [25] 
for disordered hard-core bosons in one dimension. 
After some algebra, the absorbtion rate Q becomes 



E L 



SJqItuj 



E^/ 3 [/(e a ) - f(ep)] 5(uj + e a - e ), 



where the matrix K, is defined as: 



(7) 



K aP = I E [i> a (j + + Mj)Mj + 1)] I 2 - (8) 



E„ = 



2M 



k,k< 



Vk-k' 
J 2 



[/(e*) ~ f( e k')} S(uj + e k -e k '), 

(9) 

where Vu = — ^ 2j=o e%k ^j is the Fourier transform of 
the disorder potential. Equation ^ shows that the per- 
turbation expansion in disorder is well defined provided 
the Fourier transform, Vk, is finite. 

Let us finally consider the so-called atomic limit, which 
corresponds to J <C A (yet A < [/). In this limit, tun- 
neling can be neglected and the eigenstates are given by 
"0m(i) — °~jm- From Eq. (JsJ) we find K,jji = 1 if j and j' 
are nearest neighbor and zero otherwise. The absorption 
rate, Eq. Q, thus reduces to 



L-l 



E t „ = 



'0" 

2 



-3 + 



r). 



=±1 j=0 



(10) 

In Sect. |III| we shall explicitly compare the results ob- 
tained using exact numerical diagonalization with the 
above results obtained both in the limit of weak ([9| and 
strong (10 1 disorder. Finally, it is important to empha- 



size that Eq. Q does not account for particle-hole ex- 
citations which are relevant at much higher frequencies, 
lj ~ U . These excitations become particularly important 
at unit filling (y = 1), when the strongly repulsive Bose 
gas becomes a Mott insulator and the absorption at low 
frequency u <C U predicted by Eq. ^ vanishes because 
there are no empty sites (i.e. holes) in the ground state. 
The contribution to the absorption from particle-hole ex- 
citations at unit filling will be discussed next. 
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Unit filling 



III. RESULTS FOR A DISORDER POTENTIAL 



We next turn our attention to the commensurate case 
with v = 1. For large enough U/J and U/A, the system 
becomes a bosonic Mott insulator. In Ref. it was 
shown that for clean systems the absorption rate is zero 
at low frequencies and exhibits a narrow peak of width 
~ J centered about uj = U. In this section, we consider 
the modifications of such peak due to a disorder or quasi- 
periodic potential, 

In order to obtain the energy absorption rate within 
linear response, we first use the spectral decomposition 
of the correlation function Xk{^>) i n terms of the exact 
eigenstates of the unperturbed Hamiltonian. This yields 
the following expression for the energy absorption: 



E u = 6J 2 lu^ ^\{9 n \K |* )| 2 5 (u + E - E n ), (11) 



where |\& n ) are the eigenstates of the original Hamil- 
tonian ([I]) with energies E n , and \^> ) = 1 1 , 1 , - - - , 1) is 
the ground state in the Fock representation correspond- 
ing to one boson per lattice site. The low energy states 
are particle-hole excitations \tf>(m,j)) = ^b\ n b m+ j |\&o) 
corresponding to double occupation at site m and an 
empty site m + j. These excitations are all degenerate 
with energy U in the absence of tunneling and (i.e. for 
A = J = 0). 

For finite values of A and J, the eigenstates *S> n can be 
calculated using degenerate perturbation theory by writ- 
ing = J2m,jfm,3 l^( m >i))' where the coefficients 
satisfy: 



{<t>(m,j)\H\<f>(m',j')) f m >j. = Ef„ 



(12) 



m',j' 



The above matrix element and the energy absop- 
tion can be computed by taking into account that 
(<j>(m,j) \K\^o) = — \/2<$j,±i, within the subspace con- 
taining just one particle and one hole. Hence, the ma- 
trix elements in Eq. (11) correspond to n \ K \^q) = 

~ Em,r=il V%fm,r- 

Again, let us note that, in the 'atomic limit' where 



the tunneling can be neglected, Eq. (11) simplifies con- 
siderably. Localized particle-hole excitations become 
the exact eigenstates |^ n ) = \(f>(m,j)) with energy 



E + U + e m — e m+ j, where E = J^tLi e « i s the en- 
ergy of the ground state for a given realization of . The 
absorption rate (111 then takes the form: 



M 



E u = SJqITUJ ^ ^2 <K W _ U + £ rn ~ £ m +r), (13) 



r— ±1 m— 1 



which can be readly evaluated numerically once the dis- 
order potential is known. 



In this section we assume that the on-site energy e; in 
Eq. ([I]) are random numbers uniformly distributed within 
the interval [—A, A]. The absorption rate Q can thus 
be conveniently recast as: 



E u = MSJqF, 



(14) 



where F is a dimcnsionlcss function that can be numer- 
ically evaluated. The results for incommensurate and 
commensurate cases are described below. 



Incommensurate filling 

As stated above, the absorption rate is calculated nu- 
merically at zero temperature starting from Eqs. Q and 
([8]). In Fig. [2] we plot the frequency dependence of the 
response function F for different fillings and increasing 
values of disorder A — 0.1J (upper panel), A — J (cen- 
tral panel) and A = 5 J (lower panel). Since, in the 
fermionic representation, the Hamiltonian of Eq. (pj is 
particle-hole symmetric, the absorption rate in Eq. fl7J) is 
unchanged under the transformation v — > 1 — v, so we 
restrict our discussion to fillings v ^ 1/2. 

A noticeable feature of Fig. [2] is that the response at 
high frequencies is independent of the filling factor. In 
this limit, the relevant processes contributing to the ab- 
sorption mainly involve transitions from states far below 
the Fermi level (i.e. (e a <C fi) into empty states far 
above it (i.e. [ep ^> //)). As a result, the Fermi-Dirac 
distributions in Eq. ^ become irrelevant, and thus any 
dependence on the value chemical potential fi disappears. 

However, at low frequencies, the absorption rate van- 
ishes quadratically with frequency for any strength A 
of the disorder. To understand this, let us expand 
in Eq. (Fn the distribuntion functions, f(s a ) — fisp) — 
(e a — ep)d E f(e a ). Taking into account that, at zero tem- 
perature, d E f(e) = — S(fi — e), we find that F ~ Cuj 2 , 
where the constant 



C = ^ ^ ^"/3<5(M _ £a)S(fi - Efi) 



(15) 



can be evaluated numerically. From Eq. ( 15 1 it can be 



seen that the constant C is non zero provided there are 
non- vanishing matrix elements JC a p connecting two states 
a and (3 at the Fermi level, e a = s» — fi. 

In the limit of weak disorder corresponding to A <C J, 
the absorption can be evaluated using Eq. where 
\Vk\ 2 = A 2 /3M, as follows from the expression for the 
Fourier transform of the disorder potential, Vk ■ Going to 
the thermodynamic limit and introducing the density of 
states p(e) = (2nJ v / l - e 2 /AJ 2 )-\ we find: 



F : 



uA 2 



de 



247rJ 4 y Q v /l- e 2/4J 2 v /l-(e + ^) 2 /4J 2: 



(16) 
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FIG. 2: Energy absorption rate of hardcore bosons in a ran- 
dom potential: the response function F [see Eq ( |14[ )] is plot- 
ted versus modulation frequency for different filling factors 
and increasing values of disorder strength: A = 0.1J (top 
panel), A = J and A = 5 J. Calculations are done on a ring 
of M — 500 lattice sites yielding negligeable finite size effects. 
The number of disorder realizations used was N r = 500. 



where a = max(— 2 J, p — uS) and b = min(/i, 2 J — lu). 
Here the chemical potential p is related to the filling fac- 
tor by fx = — 2</cos7tja By expanding Eq. (16) at low 
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FIG. 3: Comparison between numerics [symbols] and analyt- 
ics [Eq. Q, solid line] for weak disorder. Here A/J = 0.01 
and we consider two filling factors ^ = 0.3 and v — 0.5. The 
length of the chain is M = 2000 and the number of disorder 
realizations is N r = 2000. The inset is a zoom of the low 
frequency regime, where the absorption rate is quadratic in 
frequency as given by Eq. (171 



discussed above: 

F AV 



6J 2 



p(m) 2 - 2 



24nJ A sin 2 (7w)' 



(17) 



It should be noticed that the right hand-side of Eq. (17 1 
diverges in the limit of vanishing lattice filling v — > be- 
cause the density of states p(e) has a van Hove singular- 
ity at zero energy in one dimension. Thus the low filling 
limit, the quadratic behavior of Eq. ( 17 1 is only recovered 
at increasingly low frequencies, as shown in Fig. [2jupper 
panel) . 

In Fig. [3] we show a comparison, in the limit of weak 
disorder (A = 0.01J), of the numerical results (open sym- 
bols) obtained using exact diagonalization with the an- 
alytical expression of Eq. ( p~6[ ) (continuous lines). The 
agreement is indeed very good over the entire frequency 
range. In the inset it is demonstrated that numeri- 
cal results are consistent with the quadratic behavior of 
Eq. (17), expected at low frequencies. 



frequencies, we again obtain that the quadratic behavior 



On the other hand, in the opposite limit of strong dis- 
order, A ^S> J, the tunneling can be neglected. In this 
limit, the absorption rate can be evaluated directly from 
Eq. (10 1. Taking into account that the on-site energies 
tj at different sites are completely uncorrelated, we find 

F = 7tuj dep(e) [ de'p(e')5(u + e-e'), (18) 



where Jl and p are the disorder-averaged chemical po- 
tential and density of states, respectively. In the random 
potential the latter is constant and given by p(e) = 1/2A, 
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FIG. 4: Compari son between numerics [symbols] and analyt- 
ics [solid line, Eq |19[ )] for strong disorder. The value of the 
disorder is A/ J = 100 and the filling factors are v = 0.3 and 
v = 0.5. The length of the chain is M = 2000 and the number 
of disorder realizations is N r = 2000. In the inset we compare 



our numerics with the quadratic expansion ( 20 1 expected at 



low frequency. In both cases the agreement is quite good. 



FIG. 5: Energy absorption rate in the Mott insulator phase for 
hardcore bosons in a random potential: the response function 
F (see Eq. ( |14[ )) is plotted versus modulation frequency for 
fixed J = 0.01(7 and increasing values of the disorder strength 
A/U = 0(dashed line), 0.03, 0.06, 1. The length of the chain is 
M — 60 and the number of disorder realizations is N r = 200. 
A broadening of the absorption peak is observed for increasing 
disorder. 



and therefore fj, — (2v — 1)A. Using Eq. (18 1, the follow- 
ing low-frequency behavior is obtained: 



F : 



4A : 



: to [min(/i, A — u>) — max(— A, \x — u>)] . (19) 



This behavior is exhibited by the numerics, as shown 
in Fig. [2] (lower panel, see also discussion further below) . 
Furthermore, in the low frequency limit Eq. ( 19 ) we again 



recover the quadratic behavior described above on gen- 
eral grounds: 



F = 7X5"^ 
4A 2 



(20) 



Notice however that the proportionality constant is now 
independent of the filling factor. In Fig. Q a more de- 
tailed comparison of the numerics with the analytical re- 
sult of Eq. ( 19 1 for the limit of strong disorder is shown. 



In the inset, we also compare our numerical results with 
the quadratic behavior ( 20 I expected at low frequency. 



In both cases the agreement is remarkably good. 



Unit filling 

The absorption rate for v = 1 is calculated numerically 
starting from Eqs (11) and ( 12 1. The result is shown in 
Fig. [5] as a function of frequency for different values of 
the disorder strength A and J = 0.01E/. 

In the absence of disorder (A = 0) the absorption rate 
can be evaluated analytically and the dimensionless func- 



tion F in Eq. ( 14 ) is given by [23] : 



F 



2uj 
3J 



U 



GJ 



for \lu-U\< 6 J, 



(21) 

and vanishes otherwise. The dashed line in Fig. [5] cor- 
responds to our numerical result for A = which fully 
agrees with the above formula. 

As the strength of disorder is increased, we see from 
Fig.|5]that the absorption spectrum becomes broader and 
progressively develops a triangle-like shape as the limit 
J <C A is approached. This behavior can be obtained 
analytically starting from Eq. (13). In the J« A limit, 



since the on-site energies at different lattice sites are un- 
corrected, we obtain: 



F 



2A 2 



A nA 

de / de'5{uj-U + €-€'), (22) 

-A J -A 



where p(e) = 1/2 A is the disorder-averaged density of 
states in the atomic limit introduced above. Upon inte- 



gration, Eq. (22 I yields 



F~^(2A-\u-U\), 



(23) 



showing that the lineshape of the absorption rate in the 
atomic limit is approximately triangular and vanishes at 
\uj-U\ = 2A for J < A. 
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IV. QUASI-PERIODIC POTENTIAL 

In this section we assume that the external potential 
in Eq. is ej = Acos(2tt ja) [15j [TBI E!, being a an 
irrationalnumber. This quasi-periodic potential distri- 
bution is realized experimentally by superimposing two 
different periodic potentials with incommensurate lattice 
periods [IT] . Differently from the disorder potential con- 
sidered in the previous section, where all states are local- 
ized (in the thermodynamic limit) for an arbitrarily small 
amount of disorder, in the quasi-periodic case there is a 
phase transition [3U] at A c = 2J: for A < 2J all states 
are extended while for A > 2 J all states are exponentially 
localized. 

As described below, we find that the absorption spec- 
tra of bosons in the quasi-periodic potential are remark- 
ably different from the spectra described in Sect. |III| for 
the bosons moving on a disorder potential. 



Incommensurate filling 

For lattice fillings v < 1 (i.e. less than a boson per 
site), we calculate the absorption rate numerically start- 
ing from Eqs ^ and (JsJ) . We restrict to zero temperature 
and fix o~ — 0.77145245 which is relevant to the experi- 
ments carried out by the Florence group. 

In Fig. [6] we plot the frequency dependence of the re- 
sponse function ( 14 1 for increasing values of the disorder 
strength A. Compared to Fig. [2j we see that the absorp- 
tion spectra in quasi-periodic potential exhibit a much 
richer structure. Interestingly, the response function be- 
comes very large (and actually diverges) when u> matches 
special values. Moreover, for a given filling factor, the ab- 
sorption rate is non-zero only within a certain range of 
frequencies. 

These features can again be understood in the limits 
of weak and at strong disorder, where results can be ob- 
tained analytically. In the limit A <C J, the absorption 
rate can be calculated starting from Eq. (|9| and taking 
the thermodynamic limit. This yields (e^ = — 2Jcosfc): 

_ TTCJ [ dkdk' \V{k~k')\ 2 s . . , . 

F =^2 ^ P + + (24) 



where the integration over momenta is restricted to 
[0,27r]. Using that the (modulus square of the) 
Fourier transform of the disorder potential is |V(fc)| 2 = 
liniM^oo \ V M (k)\ 2 , where 

\v M (k)\ 2 



A^ 
4M 



1 



1 



and fcn = 2ira, for M 



gi(fc+fco) _ \ e 'i(k-k ) _ Y 

oo, the right hand-side of 



Eq. (25) becomes a series of delta functions: 

,2 +°° 



\V(k)\ = 



nA 2 
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FIG. 6: Energy absorption rate of hardcore bosons in a quasi 
periodic potential with a = 0.77145245: the response func- 
tion F [see Eq. ( 14 \] is plotted versus modulation frequency 
for different filling factors and increasing values of disorder 
strength: A = 0.1J (top panel), A = J and A = 10J. Calcu- 
lations are done on a ring of M — 1200 lattice sites yielding 
almost negligeable finite size effects. The number of disorder 
realization is N r = 500. 



Putting Eq. (26 1 this result into Eq. (24 1 we obtain: 



F 



ujttA 2 



J~] [5(k + k Q + 2im) + S(k -k Q + 2i 



dk 



4J 2 ./o 2tt 



S(oj- e k + e k +k ) 



Q(n + £ fe )6(w - efe+feo). 



(27) 



(26) The integral in Eq. ( 27 1 can be evaluated using the iden- 
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tity a sin k + bcosk = ccos(fc + 7), where c = Va 2 + b 2 
and tan 7 = —a/b. Since a = sin fco & n d & = 1 — cos ko we 
obtain c = 2sin(fc /2) and tan7 = — 1/ tan(fc /2). The 



zeroes of the 5- function in Eq. ( 27 ) occur at k = k± — 



-7 ± arccos(ix'/2 Jc). Hence, we obtain: 



F 



uj/2 



4J 2 ^(2Jc) 2 - uj< 



J2e(-s; r )e(oj + t r ), (28) 



-€k r — ^ — 2Jcosfc r — Eq. (28 1 shows 



where £ r 

that the absorption is finite only in a range of frequencies 
given by the conditions £ r < and £ r + to > 0, where 
£ r itself depens on lu through the uj dependence of k r . 
Moreover, the response diverges for uj = 2Jc, provided 
this frequency value is allowed for a given filling factor. 
Since c = 2sin7rcr = 1.31576, the divergence occurs at 
uj = 2.6315 J, as found numerically and shown in Fig. [6] 
(upper panel). 

For a potential strength comparable to the tuncling 
rate, i.e. A ~ J, the absorption spectrum develops very 
sharp peaks as shown in the central panel of Fig.[6j These 
peaks gradually disappear as A increases and becomes 
much larger than J. In this regime, however, the energy 
absorption spectrum becomes indeed rather similar to 
the case of weak quasi-periodic potential, as can seen 
by comparing the upper and lower panels of Fig. [6] This 
peculiar effect can be explained analytically starting from 
Eq. (10 1 which applies in the atomic limit. Introducing 



the variable y — 2ncrn 1 in the thermodynamic limit, we 
obtain: 



F = LUTT 



2tt 



dy 
2tt 



5(uj + A cos y — A cos(y + 2tt a)) 



e(^-Acosy)e(uj + Acosy~fi) (29) 



The integral in Eq. ( |29| ) becomes the integral in Eq. ( 27 ) 
after a change of variable k = y + it. We thus obtain: 



F 



oj/2 



6(/x+Acos k r )Q((jj— A cos k r —fi), 



r=± 



showing that the behavior of the absorption rate at 
strong quasi-periodic potential can be obtained from 
Eq. (28) by simply replacing 23 with A. 



— A=0 U 

— A=0.03U 
- A=0.06U 

A=0.1U 




FIG. 7: Energy absorption in the Mott insulator phase for 
bosons in a quasi-periodic pote ntia l with a = 0.77145245: the 
response function F [see Eq. ( 14 1] is plotted versus modu- 
lation frequency for fixed J = 0.01(7 and increasing disor- 
der strength A/U = 0(dashed line), 0.03, 0.06, 1. Here con- 
vergence is achieved for system size L — 70. As disorder 
increases, the response function broadens and changes con- 
vexity, as predicted by Eq. ( 32 1. 



does not become broad, but instead satellite peaks ap- 
pear on the sides of the central absorption feature. For 
stronger quasi-periodic potentials, the central peak dis- 
appears and the spectrum develops a two hump struc- 
ture. To understand these features, let us again focus on 
the 'atomic limit' where the absorption rate can be ob- 
tained analitically using Eq. (13 1. Introducing the vari- 
able y = 2iran and passing to the continuum limit 
obtain the result: 



we 



F = 2ujtt 



271 



dy 
2n 



S(u> — U + A cos y — A cos(y + 2jra)). 

(31) 



The integral (31 1 can be readily evaluated using the iden- 
tity a sin y + b cos y = c cos(y 
and c = "v/o 2 



& 2 . 

2 sin na and therefore, 



7), where tan 7 = —a/b 
From Eq. (31) we have that c = 



Unit filling 



We have obtained the absorption spectrum at unit fill- 



ing numerically using Eq. (Ill and Eq. (12 1, for the same 
value of a = 0.77145245. The result is plotted in Fig. [7] 
for fixed J/U = 0.01 and increasing values of the quasi- 
periodic potential strength, A. The dashed line corre- 
sponds to the clean case, A = 0, where the absorption 
rate is given by Eq. (21 1. 



We see that the shape of the absorption spectrum 
changes considerably as A increases. For sufficiently 
weak quasi-periodic potential, A < 3J, the spectrum 



F 



2u 



(32) 



showing that the absorption rate diverges at the the egde 
where \u> — U\ = Ac. This means that the shape of 
the absorption spectrum changes completely going from 
weak to strong quasi-periodic potential, as obtained nu- 
merically and shown in Fig. |jj Finally, upon comparing 
Eq. ( 23 1 and Eq. (32 1 we see that in a quasi-periodic po- 



tential the absorption spectrum is (at least in the atomic 
limit) narrower because c < 2. 
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V. EFFECTS OF A PARABOLIC TRAP 

In this section we discuss the effects of a harmonic 
trapping potential V(z) — muj 2 z 2 /2 on the absorption 
spectrum. Here m is the atom mass and a^o is the trap- 
ping frequency. In this case Vj 10 in Eq. |TJ) is non-zero and 



Quasi-periodic 



Disordered 



given by V, 



ho 



a 



ho 



(j - Af/2) 2 , where a ho = mwg d 2 /2 



d being the lattice period. 

We have repeated the calculations of the absorption 
spectrum including Vj w and the result is shown in Figj8j 
For a system of hard-core bosons in the absence of dis- 
order or quasi-periodic potential, the trap favors the for- 
mation of a Mott insulator in the center surrounded by 
a superfluid region at the trap edges. This gives rise to 
a finite absorption at low frequency (see inset in the up- 
per panel), which is related to the creation of excitations 
at the edge of the trap. By contrast, in a uniform sys- 
tem of hard-core bosons, as we have described in Sect. [II] 
the energy absorption vanishes to all orders because the 
hopping operator K commutes with the Hamiltonian. 

Let us next consider the effect of a small amount of 
disorder or a weak quasi-periodic potential. Clearly the 
Mott insulator at the center of the trap cannot absorb 
energy at low frequency, so the only contribution comes 
from the outer shell, where the filling factor is less than 
unity. In particular the low frequency peak arising from 
edge excitations fragments in multiple peaks with little 
spectral weight compared to the response from the bulk 
discussed in Sections IIIII and IIVI 

We also see in Fig. [8] that the behavior of the absorp- 
tion spectrum at frequencies close to the bandwidth cru- 
cially depends on whether the applied potential is truly 
random or quasi-periodic. Whereas the disordered case 
exhibits a smooth behavior in the absorption up to the 
bandwidth 4J where it falls to zero, the quasi-periodic 
one shows a sharp peak located at the bandwidth 4Jc 
which resembles the divergence found in the correspond- 
ing homogeneous case. Note that the position of this 
peak is almost independent on the number of atoms in 
the tube and therefore, the peak should be visible in a 
realistic experimental situation where an average over a 
multiple tube setup with variable filling is performed. 
The same conclusion applies to the system with strong 
disorder or quasi-periodic potential as can be observed 
in the lower panel in Fig. [8] For A ~ J (Fig. [8j mid- 
dle panel) the peak structure in the quasi-periodic case 
is more complex, and thus the averaging procedure will 
produce some rounding off of the peaks. Still, the absorp- 
tion can be considerably larger than in the disordered 
case and this difference should be clearly visible. 



VI. CONCLUSIONS 

In conclusion we have investigated the energy absorbed 
by a disordered strongly interacting Bose gas in the pres- 
ence of periodically modulated optical lattices. For filling 
factor less than one, the absorption rate has been calcu- 
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FIG. 8: Energy absorption rate of disordered and quasi- 
disordered hardcore bosons in a parabolic trap. The response 
function F is plotted versus modulation frequency for different 
fillings. Upper panel: A = 0.1J, middle panel: A — J, lower 
panel: A = 10J. The system size is M = 500, the number of 
disorder realizations is N r = 4000 and a h ° = 2.88 x 10~ 4 J. 
The inset in the upper panel shows the low frequency absorp- 
tion in the absence of disorder, which is finite for trapped 
gases. As A increases, this sharp peak fragments out and is 
no longer visible on the scale of the bulk contribution. 



lated exactly in the hard-core limit via the Bose-Fermi 
mapping. For commensurate filling, corresponding to 
one boson per lattice site, the gas is a Mott insulator 
and can only absorbs energy at much higher frequency 
(of the order of the repulsive interaction U). The disorder 
induced broadening of the absorption spectrum has been 
calculated by restricting to the subspace of particle-hole 
excitations. 

We have performed extensive calculations compar- 
ing two different sources of disorder: a random poten- 
tial, which is relevant for current experiments based on 
speckle patterns, and a quasiperiodic potential, which is 
obtained by superimposing two optical lattices with in- 
commensurate periods. Our results indicate that the re- 
sponse of the gas to the lattice modulation significantly 
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depends on the chosen source of randomness. 
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VII. APPENDIX 

In this Appendix we shall derive the asymptotic for- 
mula ([9]) based on perturbation theory for weak disorder. 
For clarity, we rewrite the general expression Q 



SJ§iruj 



J^/Cfcfc' [/(e*) - f(sk')] S(u + e k - e k > 



k.k' 



using new indices k and k' . Moreover we find convenient 
to introduce the matrix elements 

= E MEC? + i)VvC0 + i>W)MJ + 1)] (34) 



showing that the matrix A is diagonal in momentum 
space. Since the matrix /CjjL/ = 45^/ cos 2 k is also di- 
agonal, the absorption rate (17| vanishes. 

For A <C J, we formally expand the rhs of Eq. (34 1 
in powers of the disorder strength Ak.k' = -^-kk' 
0(A 3 ), so the matrix JC takes the form 



A 2 
A kk 



K 



k.k' 



--K, 



kk' 



\A-kk 



A-kk'i^kk' + 

,| 2 + 0(A 3 ). 



A 1 

A k'k 



A 2 

A kk< 



A-k'k) 



(36) 
(37) 



Taking Eq. (351 into account, we see that the only non- 
diagonal term appearing in the expansion (37 1 is | ^4^. fc , | 2 , 
which is second order in A. This term can be readily 



evaluated from Eq. (34) by applying first order pertur- 



bation theory for the eigenstates 



q^k 



(38) 



After a simple algebra we obtain 



A 1 - 

A kk' — 



cos k') 



(39) 



e k e k' 

which is valid up to linear order in A. Final ly, b y using 
the dispersion relation = — 2 J cos k, Eq. (39) further 
simplifies yielding 



so that Kk,k' = \ A k ,k>\ 2 - 

In the absence of disorder A = 0, the eigenstates 
are plane waves tpu(n) = e tkn /yL with energy e° = 
— 2Jcosfc. Therefore from Eq. (34 1 we find 



4° 

A kk' 



2S k k' cos k, 



(35) 



A 1 

A kk' 



J 



(40) 



Substituting Eq. ( 40 ) into Eq. ( 33 1 and replacing the 



eigenstates by their zero order values e k — e^, we recover 
the asymptotic formula 
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